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Non-magnetic ground states are a fascinating possibility allowed by the physics of 
quantum antifcrromagnets. These states - which lack the classical periodical long-range 
order - can be stabilized whenever reduced dimensionality, a small spin value, and/or 
the presence of competing interactions lead to strong enough zero-point quantum fluctu- 
ations. 

One dimensional or quasi-one dimensional spin-half Heiscnbcrg antifcrromagnets often 
have non-magnetic ground states. Indeed for the ID nearest-neighbors Heisenberg model 



the exact solution due to Bethe [1] predicts the absence of true long-range antiferromag- 
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Fig. 1. - Sketch of a spin liquid (a) and of a symmetry- broken (b) non- magnetic RVB state. 
Each stick represents a singlet bond. 

nctic order even if the ground state is very close to have a broken symmetry, with a 
gapless excitation spectrum and a power-law decay of spin-spin correlations. This is also 
the case for any array consisting of an odd number of chains (odd- leg ladder systems). 
The ground state of two chains or in general of any even-leg ladder system is non-magnetic 
too. However - in contrast to the previous cases - here the correlation length is finite 
and the spectrum is gapped. [2] Such a gap is known to decrease exponentially with the 
number of legs [3] leading to a gapless spectrum in the two dimensional limit, where the 
ground state of the Heisenberg model has genuine long-range antifcrromagnetic order. [4] 
Competing interactions may allow in principle the stabilization of a non-magnetic 
ground state even in truly two dimensional systems. One of the simplest examples of 
these frustrated systems, which has been also recently realized experimentally [5], is the 
so-called J1—J2 model [6, 7] 

(2) H = Ji Y, Si ' Sj + Ji Y S * • S j . 

n.n. n.n.n. 

where the antifcrromagnetic alignment between neighboring spins (due to J\ > 0) is 
hindered by a next-nearest-neighbors antiferromagnetic coupling (J 2 > 0). 

Classically, the minimum energy configuration of the 2D J1 — J2 model has the con- 
ventional Neel order with magnetic wave vector Q = (tt,it) for J2/J1 < 0.5. Instead 
for J2/J1 > 0.5 the minimum energy configuration is the so-called collinear state with 
the spins ferromagnctically aligned in one direction and antifcrromagnctically in the 
other, corresponding to magnetic wave vectors Q = (w, 0) or Q = (0, 7t).[8] Exactly at 
J2/J1 = 0.5 any classical state having zero total spin on each elementary square pla- 
quette is a minimum of the total energy. These states include both the Neel and the 
collinear states but also many others with no long-range order so that the occurrence of 
a non-magnetic ground state in the quantum case, for a small spin value, is likely around 
this value of the J2/J1 ratio. Indeed, at present there is a general consensus on the fact 
that the combined effect of frustration and zero-point motion leads to the disappearance 
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Fig. 2. - Variational estimate of the magnetic structure factor for the spin-half Heisenberg chain 
and two-leg ladder (filled circles). Empty dots are the numerically exact results obtained with 
the Green's function Monte Carlo method. [22] 

of the long-range antiferromagnetic order marked by the opening of a finite spin gap for 
<~ 0.4 < J2/J1 <^ 0.55. [9, 10] The nature of this non-magnetic ground state is one of the 
most interesting puzzles of the physics of frustrated spin systems. In particular an open 
question is whether the ground state of the J1 — J2 Heisenberg model is a homogeneous 
spin liquid, i.e., a state with all the symmetries of the Hamiltonian, as it was originally 
suggested by Figueirido et al. [11]. The other possibility is a ground state which is still 
SU{2) invariant, but nonetheless breaks some crystal symmetries, dimerizing in some 
special pattern (see below). [12, 13, 14, 15, 16, 10] 

A simple picture of a non-magnetic ground state can be given in terms of the so- 
called Resonating Valence Bond (RVB) states. [17] These are linear superpositions of 
valence bond states in which each spin forms a singlet bond with another spin on the 
opposite sublattice (say A and B) [18] 

(3) \iPrvb}= H r i)---H r f ) ,jf ) , 

where N is the number of sites of the lattice, r m is the distance between the spins forming 
the m th singlet bond (i m j m ), and h(r m ) is a bond weight factor. These states form in 
general a (overcomplete) basis of the S = subspace so that any singlet wave function 
can be represented in terms of them. However, they represent a non-magnetic state 
whenever the short-ranged bonds dominate the superposition (3). More precisely, it has 
been numerically shown by Liang, Doucot and Anderson[18] that the RVB state (3) has 
no long-range antiferromagnetic order for bonds that decay as rapidly as h(r) ~ r~ p , 
with p > 5. Such bonds can be either homogeneously spatially distributed on the lattice, 
with short-range correlations among each other (spin liquid) [fig. 1 (a)], or they can break 
some symmetries of the Hamiltonian, with the dimers frozen in some special patterns 
[fig. 1 (b)] as originally predicted for the J1 — J2 model in the regime of strong frustration. 
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[12, 13, 14, 15, 16, 10] 



In a seminal paper, [19] Anderson proposed that a physically transparent description 
of a RVB state can be obtained in fermionic representation by starting from a BCS-type 
pairing wave function, of the form 



where a =f, j, ek = — 2[cos/c ;r + cosk y ] is the free-electron dispersion and A& = A-k is 
the (real) gap function, provided the Fourier transform fk of the pairing function, /jj , 
satisfies: fk = Ak/(tk + + A\). The non-trivial character of this wave function 
emerges when we restrict to the subspace of fixed number of electrons (equal to the 
number of sites) and enforce Gutzwillcr projection onto the subspace with no double 
occupancies: singlet pairs do not overlap in real space and this wave function can be 
described by a superposition of valence bond states of the form (3). [19, 20, 21] 

This projected-BCS (p-BCS) wave function turns out to be an almost exact repre- 
sentation of several low-dimensional spin system with non-magnetic ground states. For 
instance it provides an excellent variational ansatz of the ground state of the Heisenberg 
chain and of the two-leg ladder giving a very accurate estimate of the ground-state en- 
ergy (fig. 2) and reproducing almost exactly the antifcrromagnetic correlations. In the 
first case the spin structure factor S(q) — (S q • S_ 9 ) shows a cusp at q = ir while for 
two- leg ladders it has a broad maximum at q = (ir,ir). These features are remarkably 
well reproduced by the (p-BCS) variational wave function (fig. 2), which generates robust 
antifcrromagnetic correlations at short distances with a very simple parameterization of 
the gap function: A fc = Ai cosfc + A 2 cos3fc for the chain and A fe = A. x cos k x -\- A.y cos ky 
for the ladder. [22] 

In two dimensions this wave function has been already studied for the pure Heisenberg 
model by several authors [20, 21] for A& oc (cosfc^ — cosk y ). In this case it provides a 
reasonable prediction for the ground-state energy but it fails in reproducing correctly the 
long-range antiferromagnetic order of the ground state. Here we show that this type of 
RVB state represents an extremely accurate variational ansatz for the J1 — J2 model in the 
non-magnetic phase when the gap function A& is carefully parameterized. In particular, 
a definite symmetry is guaranteed to the p-BCS state provided Ak transforms according 
to a one dimensional representation of the spatial symmetry group. A careful analysis 
[23, 24] shows that the odd component of the gap function A^ = — Ak + r Wt7r ) mav have 
spatial symmetries different from those of the even component Ak = &k+(-K,n)- Indeed, 
the best variational energy is obtained when the former has d x 2_ y 2 symmetry, whereas the 
latter either vanishes or it has d xy symmetry. In order to determine the best variational 





This wave function is the ground state of the well-known BCS Hamiltonian 



(5) 
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Fig. 3. - Accuracy of the ground-state energy, and overlap between the ground state and the 
p-BCS state (full dots) as a function of J2/J1, for N = 6 x 6. Empty dots are the energy 
accuracy of a Neel ordered spin-wave wave function. [26, 7] Lines are guides for the eye and the 
shaded region indicates the location of the expected transition point to the non-magnetic phase. 
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Fig. 4. - Average sign of the p-BCS state (full dots) as a function of J2/J1, for N = 6 x 6. Empty 
dots are the Marshall sign obtained with the p-BCS state with only the d x i_ y i component of 
the gap function. Lines are guides for the eye and the shaded region indicates the location of 
the expected transition point to the non-magnetic phase. 
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Fig. 5. - Antiferromagnetic structure factor of the d x 2_ y 2 (top left panel) and d x 2_ y 2 + d xy (top 
right panel) p-BCS wave functions. Lower panel: size scaling of S(ir, 7v) for the d x 2_ y 2 p-BCS 
state. 



wave function of this form we have used a recently developed quantum Monte Carlo 
technique [25] that allows to optimize a large number of variational parameters with 
modest computational effort. 

The remarkable accuracy of the p-BCS wave function in describing the ground state 
of the 2D J1 — J2 model in the regime of strong frustration can be shown by calculating 
the variational energy and the overlap with the exact ground state, \ipo), for the largest 
square cluster iV = 6 x 6 where the solution can be numerically determined by exact 
diagonalization. As shown in fig. 3 the accuracy of the p-BCS wave function rapidly 
increase by increasing the frustration ratio J2/J1 whereas conventional Neel ordered 
spin-wave wave functions [26, 7] quickly become less and less accurate. Entering the 
regime of strong frustration J2/J1 ~ 0.45 ± 0.05, where a gapped non-magnetic ground 
state is expected, the p-BCS wave function becomes impressively accurate with a relative 
accuracy on the ground-state energy of order ~ 4 x 10~ 3 and an overlap to the exact 
ground state of ~ 99%, both improved by more than an order of magnitude with respect 
to the J2 = case. This fact implies that the ground state in the strongly frustrated 
regime is almost exactly reproduced by a RVB wave function. 

Interestingly, the transition to the regime of strong frustration is marked by the 
stabilization at the variational level of a non-zero d xy component of the gap function. 
This allows to reproduce correctly the phases of the actual ground-state configurations 
as illustrated in fig. 4. A measure of the accuracy of a variational wave function \ipv) 
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Fig. 6. - Variational estimate of S(q) for (from the lower to the upper curve) N = 6 x 6, 10 x 10, 
and 14 x 14. Large empty circles are the exact diagonalization results for the N = 6 x 6. Inset: 
size-scaling of the order parameter squared. 



in reproducing the phases of the ground state can be given in terms of the average sign 
(S) = J2 X \ ( x \ t l>v)\ 2 Sg n [(x\i J v) {x\ipo)] ■ In the unfrustrated case it is well known that 
such phases are determined by the so-called Marshall sign rule [27]: on each real space 
configuration \x), the sign of the ground-state wave function is determined only by the 
number of spin down in one of the two sublattices. This feature, rigorously valid for 
J2 = 0, turns out to be a very robust property for weak frustration (J2/J1 <~ 0.3). [28] 
However, it is clearly violated when frustration plays an important role. It can be shown 
[23, 24] that the Marshall sign (i.e., (S) = 1 for J2/J1 = 0) can be obtained using the p- 
BCS wave function, with only the d x 2_ y 2 component, so that this wave function provides 
an almost exact representation of the ground-state phases for weak frustration. However, 
for J2/J1 >~ 0.4, the phases of the wave function are considerably affected by the strong 
frustration and only when a sizable d xy component is stabilized at the variational level, 
this property can be correctly reproduced. 

An even more remarkable effect associated with the d xy component of the gap func- 
tion is the change induced on antifcrromagnctic correlations. As it is shown, in fig. 5 the 
finite-size magnetic structure factor of the d x 2_ y 2 p-BCS state is sharply peaked around 
the antiferromagnetic wave vector Q = (n, tt) giving rise to a logarithmic divergence in 
the thermodynamic limit. Such a divergence is instead washed out in presence of the 
d xy component of the gap function, leading to a state with weaker short-range anti- 
ferromagnetic correlations which is of course more suitable to describe the spin-gapped 
strongly-frustrated phase. 

Of course, the accuracy in the energy does not necessarily guarantee a corresponding 
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Fig. 7. - Variational estimate of the dimer-dimer correlation functions A*'J obtained by keeping 
fixed the position of the bond (double stick) and moving the bond (k, I) (single stick) along 
the indicated patterns, d is the Manhattan distance. 6x6 (left), 10 x 10 (right). Large empty 
circles are the exact diagonalization results for the N = 6 x 6. 



accuracy in correlation functions. However, as shown in fig. 6, the comparison of the 
magnetic structure factor with the exact result gives a clear indication that correlation 
functions obtained by the variational approach are essentially exact. Furthermore using 
the stochastically implemented Lanczos technique and the variance extrapolation method 
[25] we have verified that the accuracy of the p-BCS wave function is preserved by 
increasing the lattice size. [23] 

In order to investigate the existence of long-range dimer-like correlations, as in the 
columnar or the plaquette valence bond state, we have calculated the dimer-dimer cor- 
relation functions, A i 'j = (SfSjS^Sf) — (Sf Sj)(S^,S[ ). In presence of some broken 
spatial symmetry, the latter should converge to a finite value for large distance. This is 
clearly ruled out by our results, shown in fig. 7, with a very robust confirmation of the 
liquid character of the ground state for J2/J1 ~ 0.5, which is correctly described by our 
variational approach. 

A totally symmetric spin- liquid solution proposed for this model in ref. [11] was ac- 
tually rather unexpected after the work of Read and Sachdev,[14] providing arguments 
in favor of spontaneous dimcrization. This conclusion was supported by series expansion 
[16, 10] and quantum Monte Carlo studies included the one done by two of us. [9] It is 
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clear however that it is very hard to reproduce a fully symmetric spin liquid ground state, 
with any technique, numerical or analytical, based on reference states explicitly breaking 
some lattice symmetry. 

In conclusion, the spin-liquid RVB ground state, originally proposed to explain high- 
Tcmperature superconductivity, is indeed a very robust property of strongly frustrated 
low-dimensional spin systems. Due to the success in reproducing the non-magnetic 
ground states of other low-dimensional spin systems like the ID chain and the two- 
leg ladder, [22] we expect that the p-BCS RVB wave function represents the generic 
variational state for a spin-half spin liquid, once the pairing function /jj is exhaustively 
parameterized according to the symmetries of the Hamiltonian. Work is in progress on 
this line of research. [24, 29] 

* * * 
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